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ABSTRACT 


We  describe  an  algorithm  for  minimum  delay  routing  in  a  communication 
network-  During  the  algorithm  each  node  maintains  a  list  of  paths  along 
which  it  sends  traffic  to  each  destination  together  with  a  list  of  the 
fractions  of  total  traffic  that  are  sent  along  these  paths.  At  each 
iteration  a  minimum  marginal  delay  path  to  each  destination  is  computed 
and  added  to  the  current  list  if  not  already  there.  Simultaneously 
the  corresponding  fractions  are  updated  in  a  way  that  reduces  average 
delay  per  message.  The  algorithm  is  capable  of  employing  second 
derivatives  of  link  delay  functions  thereby  providing  automatic  scaling 
with  respect  to  traffic  input  level.  It  can  be  implemented  in  both  a 
distributed  and  a  centralized  manner,  and  it  can  be  shown  to  converge 
to  an  optimal  routing  at  a  linear  rate. 
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I.  INTRODUCTION 

We  consider  the  problem  of  optimal  quasistatic  routing  of  data 
in  a  communication  network  with  slowly  varying  input  traffic  statistics. 
The  type  of  physical  situation  underlying  the  problem  is  the  one  described 
in  Gallager  [1]  and  we  refer  to  that  paper  for  further  discussion,  we 
describe  an  algorithm  for  solving  the  corresponding  optimization  problem 
in  either  a  distributed  or  a  centralized  manner.  The  algorithm  is  similar 
to  Gallager* s  method  [1]  and  its  generalized  versions  [2]  in  that  it 
relates  to  the  gradient  projection  method  for  nonlinear  programming 
(Goldstein  [3] ,  Levitin-Pol jak  [4]  -  see  Bertsekas  [2]  for  a  related 
discussion) .  The  algorithm  given  here  .is  different  primarily  in  that  it 
operates  in  the  space  of  path  flows  rather  than  in  the  space  of  link  flows. 
Furthermore  it  utilizes  a  shortest  path  computation  to  obtain  a  search 
direction  rather  than  an  upstream  summation  of  link  marginal  delays. 

As  a  result  the  amount  of  computation  per  iteration  is  typically  smaller 
for  our  algorithm.  Our  algorithm  is  also  well  suited  for  virtual  circuit 
networks  since  it  operates  directly  in  terms  of  the  quantities  of 
interest  -  the  path  flows.  In  this  respect  it  is  similar  to  a  routing 
algorithm  by  Segall  [5] .  The  algorithm  bears  also  a  close  relationship 
with  an  algorithm  developed  in  a  different  context  by  Aashtiani  [6] 
which  utilizes  second  derivatives.  Aashtiani *s  algorithm  need  not 
converge  to  an  optimal  routing  but,  for  problems  where  it  does  converge, 
the  corresponding  computational  results  have  been  very  encouraging. 

Our  algorithm  is  also  similar  to  the  flow  deviation  method  (Fratta, 
Gerla,  and  Kleinrock  [7] )  and  the  extremal  flow  method  (Cantor  and  Gerla 
[8] ) ,  in  that  at  each  iteration  it  involves  the  same  type  of  shortest 
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path  computation.  However,  the  first  of  t-hese  methods  is  characterized  by 
sublinear  and  hence  slow  convergence  rate  ( [9] ,  [10] ) ,  while  the  second  requires 
the  solution  of  a  nonlinear  programming  problem  at  each  iteration.  By  contrast 
the  convergence  rate  of  our  method  is  linear,  while  the  computational  requirements 
per  iteration  are  rather  modest.  Furthermore,  within  the  context  of  our  method 
it  is  possible  to  employ  second  derivatives  in  a  manner  that  is  reminiscent  of 
Newton's  method.  Computational  results  for  related  methods  that  employ  second 
derivatives  [11] ,  [12]  have  been  very  favorable  and  this  provides  grounds  for 
optimism  that  the  same  will  be  true  for  corresponding  versions  of  our  method. 

An  additional  benefit  of  the  use  of  second  derivatives  is  that  it  provides 
automatic  stepsize  scaling  with  respect  to  level  of  traffic  input  as  has  been 
shown  computationally  within  the  context  of  generalized  versions  of  Gallager's 
method  [11] .  This  advantage  is  of  crucial  importance  if  the  algorithm  is 
operated  in  a  distributed  manner. 

We  consider  a  network  consisting  of  N  nodes  denoted  by  1,2,...,N 
and  a  set  of  directed  links  denoted  by  L.  We  denote  by  (i,£)  the  link  from 
node  i  to  node  l  ,  and  assume  that  the  network  is  connected  in  the  sense 
that  for  any  two  nodes  m,n  there  is  a  directed  path  from  m  to  n  .  We 
consider  the  following  multicommodity  flow  problem  in  the  vector  of  total 
link  flows  f  =  {f^  |  (i,S.)e  L,  j  -  1,...,n}s 

(l)  minimize  0(f) 

subject  to  ^  f (j)  -  f  . ( j)  ■  r .  ( j) ,  V  i, j-1, . . . ,N  i  f  j 

AeO(i)  11  mel(i)  mi  1 

fu(j)  >  °,  V  (i,l).L,  i,j-l,...,N 

fjfc  (j)  ■  0  v  j  -  i,...,n  . 

N 

Here  f.,(j)  i*  the  flow  on  link  (i,£)  destined  for  node  j,  f ^  -  J]  f^tj)  is 
the  total  flow  on  link  (i,l),  0 ( i )  and  I(i)  are  the  sets  of  nodes  l  for  which 


- jJfc, 
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(i,A)£L  and  (i,i)£L  respectively/  and,  for  i  /  j,  r^(j)  is  a  given  traffic 
input  at  node  i  destined  for  j.  Each  link  (i ,1)  has  associated  with  it  a 
number  ,  referred  to  as  the  capacity  of  the  link  which  is  assumed  positive 
0*  +  «®.  The  standing  assumptions  throughout  the  paper  are: 

a)  ri(j)  >  0,  V  i, j  -  1,...,N,  i  t  j  • 

b)  The  real  valued  function  0  is  defined  on  the  set 


(2)  S  -  {f  j  0  <  fu  <  Cu  ,  V  <i,£)£L} 

and  is  convex  and  twice  continuously  differentiable.  Furthermore  lim  D(f)  =  +  00  , 

fU*CU 

—  >  0,  and  the  Hessian  matrix  of  D  is  positive  definite  for  all  fes 

3fU 

and  (i,£)eL. 

Problem  (1)  is  formulated  in  the  space  of  link  flows.  We  consider  an 
equivalent  formulation  in  the  space  of  path  flows.  Consider  the  set  of  all 


active  origin-destination  (OD)  pairs 


A  ■  {(i,j)|ri(j)  >  0,  i,j  • 

For  each  a  *  (i,j)EA  we  also  write  r. ( j)  *  r  .  For  each  OD  pair 

1  A 

«  ■  let  be  the  set  of  directed  paths  with  no  repeated  nodes 

originating  at  i  and  terminating  at  j.  For  a  “  (i,j)£A  and  a  path 

ptPa  we  denote  by  h^  the  flow  originating  at  i  destined 

for  j  and  travelling  along  the  path  p.  Let  h  *  {h  | p£P  ,  a£A}  be  the 

p  a 

vector  of  path  flows  and  H  be  the  set  of  all  feasible  path  flow  vectors 

(3)  H  *  { h I h  >  0  ,  /I*  h  »  r  V  p€P  ,  a£A  } 

P  —  JTZ  Pa  a 


Ppr  each  hEH  there  is  a  corresponding  flow  vector  f  satisfying  the  constraints 
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of  problem  (1)  and  related  to  h  linearly  by  means  of  the  equation 

(A)  f.»  -  £  £  «  (i.t)h  ,  V  (i,Jl)eL 

il  acA  P€P  p  p 

r  a 

when  5p(i,&)  -1  if  the  path  p  contains  the  link  (i,&)  and  6^ (i ,3.)  =  0 
otherwise.  Let  us  write  (4)  as 

(5)  f  -  Eh, 

where  E  is  the  appropriate  matrix,  and  let  us  consider  the  function 

(6)  J(h)  -  D(Eh) 

and  the  problem 

(7)  minimize  J(h) 
subject  to  heH  . 

It  is  easy  to  show  that  problem  (7)  is  equivalent  to  problem  (1)  in 

* 

the  sense  that  an  optimal  solution  h  of  problem  (7)  yields  an  optimal 
solution  {f ^(j) | (i,£)eL,  j  =  1,...,N}  of  (1)  via  the  relation 

f?o<j>  -  £  E  <5D(i.*)h*  * 

i  nPp  P  ” 

(i,j)eA  per(i,j) 

3D 

This  is  based  on  the  assumption  >  ®  which  precludes  the  existence 

9tU 

of  optimal  sets  of  flows  that  contain  loops. 

The  algorithm  of  this  paper  solves  numerically  problem  (7)  but  its 
definition,  as  well  as  its  convergence  and  rate  of  convergence  properties 
depend  crucially  on  assumption  (b)  made  earlier  regarding  the  function  D. 
In  the  next  section  we  describe  the  algorithm  for  the  case  where  there 
is  a  single  OD  pair.  The  generalization  to  the  case  of  several  OD  pairs 
and  the  corresponding  convergence  and  rate  of  convergence  results  are 
given  in  Section  3.  In  Section  4  we  describe  briefly  how  the  algorithm 
can  be  operated  in  a  distributed  manner. 
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Proofs  of  various  statements  and  results  of  this  paper  are  generally 
ommitted  due  to  space  limitations.  However,  the  reader  who  is  familiar 
with  recent  work  on  optimal  routing  should  be  able  to  convince  himself 
of  the  validity  of  the  algorithm  and  the  corresponding  convergence  results. 
References  where  appropriate  are  provided  as  an  aid  in  this  respect. 


2.  The  Algorithm  for  a  Single  OP  Pair 

Consider  the  special  case  of  problem  (7)  where  there  is  a  single 
OD  pair  (i,j)  *  a.  Problem  (7)  is  written  then  as 

(8)  minimize  J(h)  =  D(Eh) 

subject  to  2-i  h  =  r 

P£Pa  P  3 

*  '  •  •  •  *  «  • 

h  >  0  ,  V  ptP 

P  -  a 

The  algorithm  is  initiated  with  a  subset  of  paths  P°  C  P  ,  and  a  set 

3  3 

of  path  flows  h°  which  may  be  nonzero  only  on  the  paths  in  P°  ,  i.e., 

3 


h  =  0 

P 


if  P*P" 


and  are  such  that  the  corresponding  set  of  total  link  flows  f°  =  E  h° 
belongs  to  the  set  S  of  (2),  i.e. 


At  the  kth  iteration  we  have  a  subset  of  paths  P  C  P  and 

a  a 

k  k 

corresponding  sets  of  path  flows  h  and  total  link  flows  f  satisfying 


hk  -  0 


if  p*Pe 


fkes  . 


We  compute  a  shortest  path  p  eP  by  means  of  a  shortest  path  algorithm 

a  k 

using  as  length  of  each  link  (2.,m)eL  the  partial  derivative  . 


We  set 


pk+1  -  pku(Pk} 


4-  • 


(12) 
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(Note  here  chat  p  may  already  belong  to  Pk) .  We  then  solve  the 

a 

quadratic  programming  problem 


(13) 


minimize 


peP 


L  tdk 

k+1  p 


<V 


hk) 

P  2  p 


<V 


hk)2) 

P 


subject  to 


23  h  =  r  ,  h  >  0,  V  p£Pk+1 
0k+l  ?  3  P  " 

p£Pa 


k+1  k 

where  for  each  peP  ,  d  is  the  total  length  of  the  path  p  with 

a  p  err 


3  f  j 

respect  to  link  lengths  -r-j — L  ,  i.e. 

3flm 


(1A) 


d-  -  E  d.w,.)  , 


p  (S-.m)eL  p 


£m 


Ic 

and  y  is  a  positive  scalar.  If  h  is  the  solution  of  problem  (13) 
P  P 

we  set  for  peP 


(15) 


k+1 


if  peP 
if  p*P 


k+1 

a 

,k+l 

a 


In  the  way  of  explanation  we  mention  that  there  holds 


dk  .  k« 

p  3h  r  a 

P 

(cf.  (A),  (5),  (6)  and  (14)].  As  a  result  the  subproblem  (13)  may  be 

viewed  as  a  quadratic  approximation  of  the  original  problem  (8)  and 

the  algorithm  may  be  viewed  as  a  version  of  the  Golds tein-Levitin- 

Poljak  gradient  projection  method  (see  [  2])  -  the  only  different  being 

k+1 

that  problem  (13)  involves  only  the  paths  in  Pfl  rather  than  the 
entire  set  of  paths  Note  that  problem  (13)  involves  a  single 

equality  constraint  and  can  be  solved  by  essentially  analytical  means 
(compare  with  [2]).  It  is  possible  to  show  by  using  rather  standard 


fi'  - 
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nonlinear  programming  arguments  that  there  exists  a  scalar  y  >  0  such 
that  if 

(16)  u  <  yk  ¥  pePk+1,  k  =  0,1,... 

—  p  a 

and  the  sequences  {yk}  are  all  bounded,  then  the  generated  sequence 
k  k 

(f  }  satisfies  {f  }  CS  and  converges  to  the  unique  optimal  solution  of 
problem  (1). 

A  choice  of  y  based  on  second  derivatives  is  given  by 
P 


UP  "n 


2  k 

a  j(h  ) 

3h2 


n  S  «  <»•*> 

(m,Z)£L  BQ  p 


where  n  is  a  positive  scalar.  When  r)  =  1  the  resulting  algorithm 
represents  a  diagonal  approximation  to  Newton's  method.  Algorithms 
of  this  type  with  r)  near  unity  have  been  shown  computationally  to  be 
quite  successful  in  a  related  context  [11],  [12].  Another  possible 
choice  is 

/  *  Z  «(»•»> 

p  “  9f^  P 

k  k+1 

where  is  the  number  of  paths  in  P that  contain  link  (m,Z) . 

With  this  choice  it  may  be  shown  that  the  objective  function  of  problem 
(13)  is  an  upper  bound  to  the  one  corresponding  to  Newton's  method 
applied  to  problem  (8).  This  in  turn  can  be  used  to  show  that  if  f 
is  near  the  optimum,  then  the  kth  iteration  leads  to  guaranteed  reduction 
of  the  objective  function. 

The  algorithm  just  described  does  not  involve  a  search  for  an 

-k  k 

appropriate  stepsize  along  the  descent  direction  (h  -  h  ) ,  and  as 
such  it  is  more  suitable  for  distributed  operation.  When  the  algorithm 
is  operated  in  a  centralized  manner  it  is  possible  to  incorporate  a  linear 
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search  of  the  type  that  is  usual  in  many  nonlinear  programming  algorithms 
thereby  guaranteeing  a  reduction  of  the  objective  function  value  at 
each  iteration  (compare  with  [2]). 


The  Algorithm  for  Several  OD  Pairs 


In  the  general  case  where  there  are  several  OD  pairs  the  algorithm 
consists  of  several  simultaneous  OD  pair  iterations  of  the  type  described 
in  the  previous  section.  The  set  of  OD  pairs  A  is  partitioned  in  a 
union  of  disjoint  subsets  A^,  A2»  ...  A^  where  n  is  some  positive 
integer.  The  algorithm  operates  in  cycles  of  n  iterations.  In  the 
first  iteration  of  each  cycle  a  single  iteration  of  the  algorithm  of 
the  previous  section  is  carried  out  simultaneously  for  all  OD  pairs 
aeA^.  In  the  second  iteration  of  the  cycle  the  same  thing  is  done  for 
all  OD  pairs  in  A^,  and  so  on  until  the  OD  pairs  in  Ar  are  iterated  upon. 

At  this  point  the  n-iteration  cycle  is  completed  and  a  new  cycle  begins 
with  the  OD  pairs  in  A^. 

Among  the  possible  partitions  of  the  set  A  the  two  extreme  cases 
are  when  n  =  1,  A^  =  A,  and  when  A^  consists  of  a  single  OD  pair. 

Other  possibilities  of  interest  are  when  n  =  N  and  A^  consists  of 

all  OD  pairs  corresponding  to  the  same  origin  i  or  the  same  destination  i. 


Any  one  of  these  algorithms  can  be  shown  to  generate  sequences 
{fk}  that  converge  to  the  unique  optimal  solution  of  problem  (1)  under 
the  assumption  that  the  sequences  {y^}  are  bounded  and  satisfy  (16)  for 
a  sufficiently  large  lower  bound  y. 


The  proof  of  this  fact  is  similar  but  actually  simpler  than  Gafni's 
convergence  proof  [13]  for  a  related  algorithm  that  operates  in  the  space 


of  link  flows.  Regarding  rate  of  convergence  it  is  possible  to  show  for 


the 


case  where  n  *  1,  *  A  that  { f  }  converges  to  the  optimal  solution 


k  . 


at  a  linear  rate.  A  proof  of  this  fact  for  the  case  where  p  is  a 
constant  may  be  found  in  Bertsekas  and  Gafni  [14] .  It  is  our  conjecture 
that  the  convergence  rate  is  linear  even  when  n  >  1  but  this  remains 
to  be  shown. 


4.  Operating  the  Algorithm  in  a  Distributed  Manner 

It  is  possible  to  distribute  the  computation  involved  in  each  iteration 
of  the  algorithm  among  the  nodes  of  the  network.  Each  node  i  measures 
r.Cj)  for  all  destinations  j  and  receives  periodically  the  values  of  the 
total  flows  f^  of  all  links  (i,,m).  Node  i  can  then  execute  the 
portion  of  the  iteration  of  the  algorithm  that  corresponds  to  OD  pairs 
for  which  i  is  the  origin.  This  involves  computation  of  shortest 
paths  from  i  to  all  destinations  and  an  adjustment  of  the  flows  of 
the  currently  active  paths  according  to  (13) ,  (15) .  Each  node  i  is 
also  involved  in  monitoring  the  average  link  flows  of  outgoing  links  f^. 

The  values  of  these  flows  are  broadcast  to  all  other  nodes  at  the 
beginning  of  each  iteration  by  either  a  flooding  scheme,  or  by  means  of 
a  spanning  tree. 

A  distributed  algorithm  of  the  type  just  described  resembles  the 
current  ARPANET  routing  algorithm  [15]  in  that  information  providing  a 
length  for  each  link  is  broadcast  throughout  the  network,  each  node 
computes  shortest  paths  from  itself  to  each  destination  on  the  basis  of 
these  lengths,  and  shifts  flow  to  the  shortest  paths.  In  the  ARPANET 
algorithm  all  flow  is  shifted  to  the  shortest  path  at  each  iteration. 

In  our  algorithm  a  portion  of  the  flow  is  retained  in  previous  shortest 


paths  and,  if  this  is  done  properly,  the  resulting  flow  pattern  converges 
asymptotically  to  the  optimal.  By  contrast  an  algorithm  of  the  ARPANET 
type  cannot  provide  optimal  routing  since  it  precludes  the  possibility 
of  sending  data  along  more  than  one  paths  for  any  single  OD  pair. 

A  certain  amount  of  synchronization  of  link  flow  broadcasts  is 
necessary  in  our  algorithm.  It  would  be  interesting  to  know  what  would 
happen  if  these  broadcasts  and  the  corresponding  shortest  path  computations 
and  flow  adjustments  were  carried  out  asynchronously  as  in  the  ARPANET 
algorithm.  Such  an  algorithm  would  certainly  offer  important  practical 
advantages  but  the  associated  questions  of  convergence  are  quite  difficult 
and  as  yet  unexplored. 


-12- 


References 


[1]  Gallager,  R.G.,  "A  Minimum  Delay  Routing  Algorithm  Using  Distributed 
Computation",  IEEE  Trans,  on  Communications,  Vol.  COM-25,  1978,  pp. 73-85. 

[2]  Bertsekas,  D.P.,  "Algorithms  for  Nonlinear  Multicommodity  Network  Flow 
Problems",  International  Symposium  on  Systems  Optimization  and  Analysis, 

A  Bensoussan  and  J.L.  Lions  (eds.),  Springer-Verlag,  1979,  pp. 210-224. 

[3]  Goldstein,  A. A.,  "Convex  Programming  in  Hilbert  Space",  Bull.  Amer.  Math. 
Soc. ,  Vol.  70,  1964,  pp.  709-710. 

[4]  Levitin,  E.S.,  and  Poljak,  B.T.,  "Constrained  Minimization  Problems", 
U.S.S.R.  Comput.  Math.  Math.  Phys..  Vol.  6,  1966,  pp.  1-50. 

[5]  Segall,  A.,  "Optimal  Routing  for  Virtual  Line-Switched  Data  Networks", 

IEEE  Trans,  on  Communications,  Vol.  COM-26,  1979. 

[6]  Aashtiani,  H.Z.,"The  Multi-Model  Traffic  Assignment  Problem,"  Ph.D.  Thesis, 
Sloan  School  of  Management,  M.I.T.,  Cambridge,  Mass.,  May  1979. 

[7]  Fratta,  L. ,  Gerla,  M.  and  Kleinrock,  L. ,  "The  Flow  Deviation  Method: 

An  Approach  to  Store-and- Forward  Communication  Network  Design",  Networks. 
Vol.  3,  1973,  pp.  97-133. 

[8]  Cantor,  D.G. ,  and  Gerla,  M. ,  "Optimal  Routing  in  a  Packet  Switched 
Computer  Network",  IEEE  Trans,  on  Computers,  Vol.  C-23,  1974,  pp. 1062- 1069. 

[9]  Cannon,  M.D.  and  Cullum,  C.D.,  "A  Tight  Upper  Bound  on  the  Rate  of 
Convergence  of  the  Frank-Wolfe  Algorithm",  SIAM  J.  Control,  Vol.  6, 

1968,  pp.  509-516. 

[10]  Dunn,  J.C.,  "Rates  of  Convergence  of  Conditional  Gradient  Algorithms 

Near  Singular  and  Nonsingular  Extremals" ,  SIAM  J.  Control  and  Optimization 
Vol.  17,  1979,  pp.  187-211. 

[11]  Bertsekas,  D.P.,  Gafni,  E.  and  Vastola,  K.S.,  "Validation  of  Algorithms 
for  Optimal  Routing  of  Flow  in  Networks",  Proc.  of  IEEE  Conf.  on 
Decision  and  Control,  San  Diego,  Calif.,  Jan.  1979,  pp. 220-227. 

[12]  Vastola,  K.S.,  "A  Numerical  Study  of  Two  Measures  of  Delay  for  Network 
Routing",  M.S.  Thesis,  Department  of  Electrical  Engineering,  Univ.  of 
Illinois,  Urbana,  Ill.,  Sept.  1979. 

[13]  Gafni,  E. ,  "Convergence  of  a  Routing  Algorithm",  Lab.  for  Information 
and  Decision  Systems  Report  R-907,  M.I.T.,  Cambridge,  Mass.,  May  1979. 

[14]  Bertsekas,  D.P.  and  Gafni,  E.,  "Projection  Methods  for  Variational 
Inequalities  with  Application  to  the  Traffic  Assignment  Problem", 

Lab.  for  Information  and  Decision  Systems,  Working  Paper,  M.I.T.,  Cambridge, 
Mass.,  Feb.  1980. 

[15]  McQuillan,  J.M. ,  Richer,  I.,  and  Rosen  E.C.,  "The  New  Routing  Algorithm 
for  the  ARPANET",  Bolt,  Beranek  and  Newman,  Inc.,  Report,  May  1979. 


Distribution  List 


Defense  Documentation  Center  12  Copies 

Cameron  Station 
Alexandria,  Virginia  22314 

Assistant  Chief  for  Technology  1  Copy 

Office  of  Naval  Research,  Code  200 
Arlington,  Virginia  22217 

Office  of  Naval  Research  2  Copies 

Information  Systems  Program 
Code  437 

Arlington,  Virginia  22217 


Office  of  Naval  Research  1  Copy 

Branch  Office,  Boston 

495  Sumner  Street 

Boston,  Massachusetts  02210 

Office  of  Naval  Research 
Branch  Office,  Chicago 
536  South  Clark  Street 
Chicago,  Illinois  60605 

Office  of  Naval  Research 
Branch  Office,  Pasadena 
1030  East  Greet  Street 
Pasadena,  California  91106 


Naval  Research  Laboratory 
Technical  Information  Division,  Code  2627 
Washington,  D.C.  20375 

Dr.  A.  L.  Slafkosky  1  Copy 

Scientific  Advisor 

Commandant  of  the  Marine  Corps  (Code  RD-1) 

Washington,  D.C.  20380 


1  Copy 


1  Copy 


6  Copies 


Office  of  Naval  Research 
Code  455 

Arlington,  Virginia  22217 


1  Copy 


Office  of  Naval  Research  1  Copy 

Code  458 

Arlington,  Virginia  22217 

Naval  Electronics  Laboratory  Center  1  Copy 

Advanced  Software  Technology  Division 
Code  5200 

San  Diego,  California  92152 

Mr.  E.  H.  Gleissner  1  Copy 

Naval  Ship  Research  &  Development  Center 
Computation  and  Mathematics  Department 
Bethesda,  Maryland  20084 

Captain  Grace  M.  Hopper  1  Copy 

NAICOM/MIS  Planning  Branch  (0P-916D) 

Office  of  Chief  of  Naval  Operations 
Washington,  D.C.  20350 


Advanced  Research  Projects  Agency  1  Copy 

Information  Processing  Techniques 
1400  Wilson  Boulevard 
Arlington,  Virginia  22209 

Dr.  Stuart  L.  Brodsky  1  Copy 

Office  of  Naval  Research 
Code  432 

Arlington,  Virginia  22217 

Captain  Richard  L.  Martin,  USN 
Commanding  Officer 
USS  Francis  Marion  (LPA-249) 

FPO  New  York  09501 


1  Copy 


